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Abstract 

We present exact results for the dynamical structure function, i.e. the density- 
density correlations for the 1/r 2 system of interacting particles at three special 
values of the coupling constant. The results are interpreted in terms of exact 
excitations of the model which are available from Bethe's Ansatz, thereby 
throwing light on the quasi-particle content of the elementary excitations. 
We also obtain the first moment of the discrete version of the model, which 
provides a non-trivial check on its structure function. We show that the 

property of spectrum saturation is a common feature of both versions. 
03.65.-w,75.10.Jm 
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I. INTRODUCTION 



The 1/r 2 system of interacting particles, introduced by Calogero |T[ and Sutherland 0, 
and indeed even earlier by Dyson within a relaxational framework in the course of his 
discussion of the Brownian motion of Random Matrices, has continued to be of great interest. 
In the recent past, it has also generated great interest in the context of a discrete version, i.e. 
a spin 1/2 model introduced independently by Haldane and Shastry KJ, and in terms of its 
algebraic content PHT(]|j. Its interest derives from the combination of beautiful mathematical 
structure and rich physical phenomena of Quantum Fluctuations in a low dimensional system 
(quasi-LRO, non Fermi Liquid behavior etc.), as well as surprising tractability. 

In a recent development in an apparently completely different physical system, namely 
that of electrons in a random medium, Simons et al. [0,0 have succeeded in computing a 
certain correlation function depending on two variables, say space and time, and conjectured 
that this represents the density- density correlation function of the above 1/r 2 model at three 
appropriate values of the coupling constant (3 = 1, 2, and 4, corresponding to orthogonal, 
unitary, and symplectic ensembles, respectively. The results are obtained due to a suggested 
equivalence of the problem to the evolution of energy eigenvalues of a disordered metallic 
grain subject to an arbitrary perturbation |13[ to that of the 1/r 2 many body problem. The 
mapping, performed on the level of the two-point, time-dependent, density-density correla- 
tion functions, leads to an explicit exact result for the density-density correlation of the 1/r 2 
model for the above values of the coupling constant. The result, astonishingly enough, is 
valid at all length and energy scales, not just in asymptotic regions. This conjecture has been 
explicitly confirmed in a recent work by Narayan and Shastry |f[4 , where they established 
the correspondence between the evolution of the distribution of eigenvalues of a random 
matrix subject to a random Gaussian perturbation, and a Fokker-Planck equation which 
is equivalent to the 1/r 2 model. At the same time, Simons et al. [H| have established a 
direct connection between the 1/r 2 model and correlations in the spectra of random matrices 
through a continuous matrix model. 
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The excitation spectrum of the 1/r 2 model is available in great detail from the (Asymp- 
totic) Bethe's Ansatz (ABA) invented by Sutherland |IB| . The picture that arises is that 
of an underlying gas consisting of quasi-particles obeying Fermi statistics, and interacting 
weakly with each other through a Hartree-Fock interaction leading to a back flow. This 
picture indeed gives all the excited states of the model. The remaining problem then, is 
that of an appropriate decomposition of the "bare" particles into "quasi-particles". The ex- 
plicit knowledge of the correlation functions provides us with an opportunity to describe the 
intermediate states in S(q, uS) phenomenologically as combinations of the "quasi-particles" , 
whose energies are available from the ABA. This is analogous to quark spectroscopy in the 
theory of elementary particles, the ABA quasi-particles and quasi-holes are our quarks in 
the present scheme. We also provide a nontrivial check on the conjecture by comparing the 
explicitly known functions of the momentum q of the 1/r 2 model, with those obtained by 
integrating the explicit correlation functions. 

The plan of the paper is as follows. In Section [TI|, we define the 1/r 2 model, and sum- 
marize the known information about its moments. The dynamical correlations found in 
Refs. [0,0] are summarized and their Fourier transforms given. The small q Hydrody- 
namic limit of these correlators is calculated, and the saturation of the spectrum by a sound 



like linear mode (i.e. 00 = sq) is demonstrated. Section |T| contains a summary of the results 
of the Asymptotic Bethe's Ansatz for the 1/r 2 model, where we write down the dispersion of 



the effective quasi-particles and quasi-holes. In Section |TV|, we rework the expressions for the 
structure function for the three ensembles into forms wherein the energy conserving delta 
functions are shown to have a natural interpretation in terms of multi quasi particle-hole 
pairs. In Section [V|, we discuss the discrete 1/r 2 model, and display its first three relevant 
moments explicitly. The similarity to the continuum 1/r 2 model in terms of the exhaustion 
of the structure function by "spinons" at low q is pointed out, and this is highlighted to 
be a unique feature characterizing this family of models. In Section |VT], we summarize our 
results. 
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II. UNIVERSAL CORRELATION FUNCTIONS APPLIED TO THE 1/R 2 MODEL 



The Sutherland-Calogero-Moser system with periodic boundary conditions has a Hamil- 
tonian describing spinless fermions confined to a ring and interacting through a 1/r 2 pairwise 
potential [0: 

For simplicity we have chosen the mass to be 1/2. The ring has length L and the number 
of particles is N. The statistics of the particles can be chosen arbitrarily since the particles 
cannot get past each other owing to the singular nature of the interaction at the origin; 
we declare them to be fermions for convenience (they are indeed so at /3 = 2), and one 
might imagine the system to be that of fermions with either repulsive {(3 > 2) or attractive 
((3 < 2) interactions. Unless explicitly specified, we will assume that the system is in the 
thermodynamic limit (L — > oo and iV — > oo) and has a finite 0(1) density (d = N/L). 
It was first argued in Ref. |ll| that the time-dependent correlation functions of this 



one-dimensional Hamiltonian are equivalent to certain universal correlation functions JT3 
of the energy spectra of weakly disordered metallic grains when (5 = 1, 2, or 4. For spectra 
subject to some arbitrary perturbation, X, exact analytical expressions were derived for the 
two-point density correlation function, 

k(E, X) = (p(E -E,X + X)p(E, X)) - (p(E, X)) 2 , (2) 

where p(E,X) = J2{=i 5(E — Ei(X)) is the density of states of the system and (• • •) denotes 
a statistical average which can be performed over a range of energy or over X. It was shown 
that after the following rescaling in which the parameters become dimensionless, 6j = Ei/A 



and x = X y {(dei(X) / dX) 2 ) , Eq. (g) becomes universal, depending only on the symmetry 
of the Dyson ensemble. In fact, the universality is not specific to disordered metals but 
applies equally to all non-integrable or quantum chaotic systems |fT~3f| . 
Remarkably, by performing the change of variables 
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x 2 = -2it ,e = r , (3) 

where t is the time coordinate and r is the spatial coordinate, Eq. (0) becomes equivalent 
to the two-point particle density correlator of the ground state of the Sutherland model. 
The coordinate r will be given in units of the mean interparticle distance 1/d (although 
d = (p) = 1, we will continue to display the "d" dependence in order to retain generality). 
The resulting correlation function, after the change of variables, is 

k(r,t) = (p{f-r,i+t)p{f,i))-d 2 , (4) 

with p(r,t) = Y,iLi <K r — r i{t))- The expectation value (• • •) is to be taken on the ground 
state of H. 

Once we have an expression for k(r,t), we can calculate the dynamical structure factor 
S(q, ui) |T7|. The explicit connection between k(r, t) and S(q, ui) is made by taking the space 
and time Fourier transforms of k(r,t), 

S(q,u) = — fdr [dt k(r,t) e^'^ . 



2nd. 

S(q, uj) has a representation in terms of the excited states of the system: 

= jjT, KHp,|o>I 2 ^ -e v + e ) , (5) 



where H\v) = E u \u), and 



r N 
p g = dr p(r) e~ iqr = ^ e~ iqn . (6) 



In the following, we will present exact analytical expressions and discuss some of the 
properties of S(q, oj) for the three special values of j3. 

A. The Moments of S(q,uj) 



We begin by stating some important sum rules [17 1 concerning the function S(q,u). 
Defining the moments of this function as 
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In(q) = du u n S(q,u) 
Jo 



it follows from the velocity independence of the interaction that 



(7) 



h{q) = q 2 



(8) 



This is the statement of particle conservation and is the familiar /-sum rule. Another sum 



rule follows from the compressibility relation [17 



lim I-i(q) 



q^O 



(9) 



where s 2 = 2(dP/dd) is the square of the sound velocity and P is the pressure. However, 
since P = —(dE/dL) and for the Sutherland model the ground state energy is known to be 



E = (tc 2 f3 2 / 12) (N 3 /L 2 ) [gj, it follows that Eq. <g) can also be written as 

lim J_i (g) = ) . 

q^O ^ Z d Z p Z 



(10) 



Finally, we note that the zeroth moment, Io(q), also called the static form factor (usually 
denoted by S(q)), has been determined before for all three values of (3 from Random Matrix 
Theory (RMT) [)T8| . The connection between the 1/r 2 model in the static limit and the 
distribution of eigenvalues of a random matrix was established by Sutherland in his early 
papers ||. Therefore, from the results of Ref. [§] we anticipate the following expressions: 
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io(q) 



2k, 



-[\q\ + (2k F -\q\)e(\q\-2k F )\ {(3 = 2) 



(11) 



and 



I°o(q) 



M 

k p 



l->(l + g) 

\q\+k F 



\q\-k F 



io(q) 



ik F 



| In 



1 _ _M_ 

2k F 



For small q we find the limiting behavior 



\q\ < 2k F 
\q\ > 2k F 

\q\ < ^k F 

\q\ > Ak F 



Io(q) 



q ^o _q_ 



09 = 1) 



09 = 4) . 



(12) 



(13) 



(14) 



In terms of the three moments given here, we can calculate two characteristic frequencies 
that give us an idea of the dispersion relations of the excited modes, namely the Feynman 



spectrum pjl, 



^ (5) = W) = W) ' (15) 

which has the small q behavior 

Mq) *3 Pfa , (16) 

and another one we call the "hydrodynamical spectrum" , 

Mi) = k ^ fak F . (17) 

In Fig. 1 we have plotted these dispersion relations for the three values of (3 for which we 
know the moments exactly. For the repulsive and noninteracting cases there is a logarithmic 
dip at q = 2kp, while for the attractive case the dispersion grows monotonically. The 
appearance of a dip can be interpreted as a tendency towards "crystallization" (i.e. the 
particles tend to arrange themselves in a lattice with spacing 1/d) as the strength of the 
repulsive interaction increases. 



B. Static Correlation Functions: Real Space 

We note that the static density-density correlations are simply related to the Fourier 
transforms of the moments Io(q). Writing the density correlation function in the form 

k(r,0) = d5{r) + d 2 C{r) , (18) 

the dimensionless correlation function C(r) satisfies the relations lim r ^oC( r ) — > — 1 and 
lim r _ >00 C(r) — > 0. It may be written as 

C7(r) = i/_" ^oos(gr)[/ (ff)-l] . (19) 

The correlation function, in scaled variables, has the representation 
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C(f/d) = r -1 cos(-Kqr)[I Q (qk F ) - 1] . (20) 

J -co Z 

Explicit expressions are available for the various ensembles from Mehta |20| by noting that 
C{r = r/d) is nothing but the two- level cluster function — ^(r), where f is the separation 
in units of the average interparticle spacing l/d. 



C. Dynamical Correlations 



We now recapitulate the results from Ref. |n|]i"2|j for the dynamical correlation function 



and present explicit expressions for S(q,u). We note that S(q,u) is real and positive; 
moreover, it vanishes for u < and it depends only on the absolute value of q. 



Unitary Ensemble 

We will first examine the simplest case when (5 = 2 and the system is non-interacting, 
which corresponds to the unitary ensemble. It can be readily shown that 

k u (r, t) = —J d\ x J i d\ exp[-ik 2 F t{Xl - A 2 )] cos[£; F r(Ai - A)] , (21) 

where A = k/k F , Xi = ki/kp, and kp = nd (kp is the Fermi momentum). Taking the Fourier 
transform in both space and time we get 

1 rco rl 

S u (q, u) = 2tfl dX i < /_ 1 dX - A ' " " A " 

' 6{u + q 2 - 2k F \q\) 6(u - q 2 + 2k F \q\) 8{2k F \q\ + q 2 - u) . (22) 



Ak F \q 



In the Fig. 2 we have plotted the region of support corresponding to Eq. (0), which is 
nothing but the particle-hole continuum (in this case all excitations are in the form of 
pairs). The tridimensional plot of S u (q,u) is shown in Fig. 3. 

From Eq. (p2"|) we can of course compute all the moments of S u (q,u) exactly. The three 
moments J", Jg , and J" 1 are plotted in Fig. 4 and we remark that they are in agreement 
with the sum rules of Eqs. (H,0) and the identity Eq. fllTl). 
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Orthogonal Ensemble 



Secondly, we will consider the orthogonal (attractive) case, when (3=1. This value of 
the coupling constant leads to a more complicated expression for the two-point correlation 



function; after Ref. [[L1J| , we have 

x exp[-ik 2 F t(2\l\ 2 2 - A 2 - A 2 , - A 2 + l)/2] cos[A; F r(A 1 A 2 - A)] . (23) 

Taking the Fourier transform of k°(r,t) in space and time yields 

co, n V/- 00 ,, r» [1 - (AiAg - |g|/M 2 ] 

S ^ Uj) -^k dXl Ji dX2 (A? + A] + g 2 /4 - AfAl - l) 2 

x 5(A 2 + A 2 + g 2 /^ - 1 - A 2 A 2 - 2AiA 2 |g|/A; F + 2w/Jfe|.) 

x ^(AiA 2 - \q\/k F + 1) 0(1 - AiA 2 + \q\/k F ) . (24) 

The region of support of S°(q, uj) is plotted in Fig. 5. and it shows a continuum similar to the 
excitation of a single particle-hole pair in the non-interacting case. However, as we will later 
discuss in detail, the excited states have a more complex structure. The tridimensional plot 
of S°(q, uj) obtained by numerically integrating Eq. (^4j) and using a Gaussian regularization 
of the delta function is shown in Fig. 6. The ridge along u = q{q + kp) indicates an algebraic 
(inverse square root) divergence. 



The evaluation of Jq starting from Eq. fl24|) has been done analytically by Efetov |2T 
and it can be readily generalized for any 1°, n > 0. However, we have only been able to 
evaluate numerically the moments with n < 0. The moments 7°, Jq, and I°_ x are plotted in 
Fig. 7 and they obey exactly the sum rules of Eqs. (||) and (|K]), and the identity Eq. (|12|). 

Symplectic Ensemble 

Finally, we look at the repulsive case, when (3 = 4. In the context of RMT this cor- 
responds to the symplectic ensemble. We start with the correlation function originally 
obtained in Ref. fOl, 



kHr t) - ^ (°°dX f 1 dX f 1 dX (A 2 - 1)(A - AjA 2 ) 2 

x exp[-4i£|t(A 2 + X\ + A 2 - 2A 2 A^ - 1)] cos[2/c F r(A - AiA 2 )] . (25) 

We take the Fourier transform of k s (r, t) to get 

qs{ s <f /* f 1 » [(AiA 2 + |g|/2fc F ) 2 -l] 
5 (g ' w) ~ 64fc£ i-i dXl L A2 (A 2 + A 2 + g y4P F - A?A 2 - 1)2 

x 5(Xl + A 2 + g 2 /44 - 1 - A 2 A 2 + AiA 2 |g|/A; F - uo/Ak 2 F ) 

x 0(AiA 2 + \q\/2k F - 1) . (26) 

The region of support of S s (q,u) is shown in Fig. 8. The continuum reaches u; = not only 
at q = and 2/cp, but at q = 4kp as well. The figure is similar to the one obtained when we 
excite two particle-hole pairs in the non-interacting case; in fact, the real structure of the 
excited states is more complicated than this simple picture, as we will later demonstrate. 
In Fig. 9 we show the tridimensional plot of S s (q, uj) with a Gaussian regularization and 
a numerical integration of Eq. (p6|) . The ridges along uj = 2q(2kp — q) and uj = q 2 — 4k F 
indicate an algebraic (inverse square root) divergence. 

As for the orthogonal ensemble, the analytical evaluation of J* can be done for n > by 



generalizing the method of Ref. [^T|] for n — 0. For n < we have only been able to proceed 



with a numerical evaluation. In Fig. 10 we have plotted the moments If, lg, and Jf^. It is 



simple to check that the moments obtained in this way also agree with Eqs. (|8|,p!0|JT3|) . 



D. Spectrum Saturation in the Hydrodynamic Limit 

An important property of the Sutherland model is the saturation at q — > 0. As we have 
pointed out before in Eqs. ( P^ , |l7|) , both the Feynman and hydrodynamical spectra tend to 
the same value in this limit. In fact, this is also true for any characteristic frequency defined 
as the ratio of any two moments. We shall prove it for the orthogonal ensemble: In Eq. (p4|) 
we change variables to 

Ai )2 = 1 + ^1,2, 
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which yields 



n+l 



2q 2 ( k'jA 



kp 



■» [1 - (l + x+ +x lX2 -q/k F f\ 

(lX\ I &X2 7 — / 9 TTTT 



x 



k f 



+ Xix 2 (A + x x x 2 + 2x + ) 



x 9(2 + x + + x\x 2 — q/k F ) 9(q/kF — x + — x±x 2 ) 



(27) 



where x± = X\ ± x 2 . This last expression can be much simplified in the limit q — » 0; it 
becomes 



fcp JO jO 



rq/k F -x\ 



(q/k F - z+) 



(g 2 /4-4 + x 2 ) 2 • 

Changing the integration variables to x± and performing the double integral we obtain 



(28) 



\n+l 



(29) 



As a result, 



*n+m(g ) 



1 1/m 



qk F . 



(30) 



The proof for the symplectic and unitary ensembles is quite analogous; one obtains for the 
three values of (3 



in(q) -^(/?gM n+1 



(3 2 k 2 r 



(31) 



This yields, for any integers n and to, 



In+m(q) 

XM) 



1/m 



(3qk F . 



(32) 



The saturation property can also be visualized in Figs. 2, 5, and 8: the fact that the 
lower and upper branches of parabola have the same linear term as q — > implies that all 
characteristic frequencies must have the same asymptotics. 
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III. ASYMPTOTIC BETHE'S ANSATZ 



We now summarize the results of the asymptotic Bethe's Ansatz, which gives an explicit 
expression for the "particle-hole" like excitations underlying the system. The complete 
excitation spectrum for the 1/r 2 model can be described in remarkably simple terms as 
follows: The total energy of a state of the system is expressible as 

E = £ Vl , (33) 

n 

with the "pseudo-momenta" p n satisfying the equation, 

Pn = k n + —4^ ^ £ si S n {k n ~ k m ) . (34) 
The total momentum of the state is 

P = £ = £ • (35) 

n n 

The bare momenta are given by k n = 2nJ n / L, where the J n 's are fermionic quantum numbers 
Ji < J 2 < J3 ■ ■ ■ < Jn- Note that at (5 = 2 the interaction is turned off and we recover the 
free-fermion results. The important point is that the totality of states for the N particle 
sector is obtained by allowing the integers J n to take on all values consistent with Fermi 
statistics, not only for ft = 2, but for all ft G [1, +00]. The summation in Eq. ( |3~3| ) is trivial 
to carry out and we find 

(ft - 2)7T / N + l\ 

Pn = k n -\ n — . (36) 



L 

We can now select an arbitrary state of the system by specifying that states {ki, k 2 , . . .} are 



occupied, i.e. by introducing the fermionic occupation numbers n(kj) = 0, 1, such that 
E = ^e(k n )n(k n ) + v (kn ~ k m )n(k n )n(k m ) + l> 



(37) 



with e(k) = k 2 and v(k) = ir(ft — 2)\k\/2L. For future reference, the ground state is 
represented by no(k n ) = 1 for \k n \ < kp and n Q (k n ) = otherwise, where Uf = nd. We 
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remark that the above expression of the energy takes the form of a renormalized Hartree- 
Fock theory; a Hartree-Fock energy expectation value of the interacting Hamiltonian in a 
determinantal state n C I-|0) leads to precisely this type of expression. Note that the Fourier 
transform of the two-body interaction can be deduced from the expansion 

Therefore, Eq. ([57]) states that a Hartree-Fock expression with a renormalization of the 
coupling constant /3(/3 — 2) — ► 2(j3 — 2) leads to the exact spectrum of the model. 

We now consider the excitation spectrum near the ground state, wherein we excite a 
particle-hole pair in the free Fermi system and ask what the energy of the interacting system 
is by including the Hartree-Fock back flow term. From this point onwards we measure all 
momenta in units of hp and energies in units of the Fermi energy. Let us suppose that one 
of the particles described by Eq. (|3~?D has initially a momentum k, with \k\ < 1; we promote 
it to some state labeled by k + q, with \k + q\ > 1. The energy cost in units of the Fermi 
energy is equal to 

A(q, k) = e{k + q)- e(k) + 2 ]T [v(k + q - k') - v(k - k')} 

|fe'|<i 

= q 2 + 2kq + ^^-(2\k + q\ - k 2 - 1) , (38) 

and the momentum of this state is simply q. This implies that we can associate a generalized 
energy corresponding to a particle e > (k) (i.e. \k\ > 1) and a hole £<(/c) (i.e. \k\ < 1): 

e>(k) = k 2 + (/3 - 2)\k\ 

^<W = ^ 2 +(f-l) , (39) 

such that 

A(q,k) =e > (k + q) -e<(k) . (40) 

Note that e > (k) and e < (k) are continuous and have continuous derivatives across the Fermi 
surface. These expressions for the particle and hole energies look different from the results of 
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Sutherland ]16[ but their equivalence may be readily checked. We have chosen to label our 
states by the "bare" momenta u k n " , while Ref. [|16j works with the pseudomomenta "p„"; of 
course these are in one to one correspondence and so the choice is a matter of convenience. 

We will introduce in the usual way, particle operators A\k) and hole operators B\k) 
with the convention that the momenta corresponding to these are constrained by \k\ > kp 
for particles and \k\ < kp for holes, with excitation energies 

E A (k) = £>(&) - fi 

E B (k)= f i-e < (-k) , (41) 

where p = £>(/cf) is the "chemical potential". The quasi-particle quasi-hole excitation 
created by the operator A^(k + q)B^(-k) then has energy E A {k + q) + Es{—k), which of 
course is equal to A(k, q). Having introduced the underlying fermionic quasi-particles quasi- 
holes through Eqs. (|39|), we would like to see if the excitations generated by the bare density 
fluctuation operator p q can be expressed in terms of the latter. One of our objectives then, 
is to express the excitations of the system probed by the bare density fluctuation operator 
p q in terms of the quasi-particle quasi-hole operators. Recall that in Landau's Fermi Liquid 
Theory P3 one expresses the bare particles c(k) in a series involving quasi-particles and 
quasi-holes of the form 

c(k) = y/T k B\-k) + ]T M[k,p, 1} B\p)B\l)A\-k - p -/) + ... , (42) 

(p,0 

where \k\ < kp, and a similar expansion for particles, where z k is the quasi-particle residue. 
The density fluctuation operator p q = Y^,k c Kk + q)c(k) then has a development in terms of 
1,2,3,... pairs of (quasi) particle-hole excitations. In one dimension, we expect z k to vanish 
for arbitrary non-zero interactions, and hence the particle-hole series is expected to be such 
that the single pair should not appear. The expansions are somewhat non-unique, in view 
of the fact that we can add an arbitrary number of "zero energy" and "zero-momentum" 
particle-hole excitations to any given scheme. 
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IV. QUASI-PARTICLE CONTENT OF THE STRUCTURE FUNCTION 



For the unitary case there is no interaction and consequently quasi-particles and quasi- 
holes are the same as particles and holes: Eq. (^) at (3 = 2 exactly describes the spectrum 
of Fig. 2. On the other hand, for the orthogonal and symplectic cases the simple creation 
of quasi particle-hole pairs cannot account for the whole excitation spectrum. In order to 
see that in general (/3 ^ 2), we begin by considering the spectrum (q x u) for a particle- hole 
excitation (Eq. fl38|)), which is the familiar pair spectrum renormalized by the interaction. 
For a fixed q, the maximum value of u; occurs when k — 1: io max = q 2 + f3\q\. The minimum 
value of u depends on q: for |g| < 2 it occurs when k = 1 — \q\; for \q\ > 2 it occurs when 
k = -1. This results in u min = /%|(2 - \q\)/2 for q < 2, and u min = (\q\ - 2)[\q\ + ((3 - 2)], 
for q > 2. The curves bound a continuum which does not agree with the hashed regions of 
either Figs. 5 or 8. 

We can also promote two quasi-particles from the Fermi sea, instead of just one. The 
result is that the upper limit for uj is then given by + (3) and the lower limits are 

/%|(2- \q\)/2, /3(|g| -2)(4- |g|)/2, and [\q\ - 4)[|g| + 2(/3 - 2)]/2 (for the intervals \q\ < 2, 
2 < \q\ < 4, and \q\ > 4, respectively). Again, we note that the continuum bound by these 
curves is not equal to the hashed region in Fig. 8. 

Below we recast the expressions for the correlation functions in the various ensembles in 
terms of new variables, in order to reveal their exact quasi-particle content. 



Orthogonal Ensemble: Change of Variables 



We now turn to the expression Eq. ([24] ) of the structure function. Firstly we change 
variables and introduce 

u = AiA 2 

z = Ai + A 2 . (43) 
With this change of variables, we find 
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S°(q, u) = 2q 2 / du / dz == [l-(u- q) 2 ]5(E)/D , (44) 

Jmax{l,q-1} J2yJH yZ z — Au 

where 

s/D = z 2 + q 2 - (1 + u) 2 
E = \[~D + 2to - 2uq . (45) 

In the expressions above and hereafter in this section we will set q > without loss of 
generality. We change the integration variable u by defining k = u — q, in terms of which 
the k integration is restricted to max{l — q, — 1} < k < 1, and is immediately recognizable 
as the momentum of a hole restricted to the Fermi surface with \k + q\ restricted to be a 
particle. The z integration can be conveniently transformed by introducing a "rapidity" 
variable 

z = 2^(k + q) coshfl , (46) 

and recalling that for the orthogonal case £>{k) = k 2 — \k\ and e < (k) = (k 2 — l)/2 (see 
Eq.©), so that 

oo/ x <f f j h [~g<(fc)] ^v 7 ^) 
2 Jfc+ g >i [o; - q{k + q)\ 2 Jo 

x 5(e>(A; + g) - £<(&) - 2(jfe + q) sinh 2 - u) . (47) 

The 6 1 integration can be done simply, and gives the final result 

q 2 r [- e< (AQ] go; - A(fc, g) + (fc + g - l) 2 /2) fl(A(fc, g) - u) 
S (q,u) = — | fe |<i dk — — -f- = — = . 

4 Wi [^-<?(^ + <?)] 2 ^/A{k,q)-u ^A(k,q) + 2(k + q)-u 

(48) 



Bethe Quasi Particle-Hole Content: Orthogonal Ensemble 

We can rewrite the energy conserving delta function in Eq. ([17]) as 5(e>^(A; + g|#) — 
£<(&) — a;), where e>\k\9) = £>(k) — 2A;sinh 2 6 l . One possible picture suggested then has 
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the excited state particle possessing a "hidden" gauge variable 9, which lies in a limited 
range, as a hole would, endowed with energy but possessing no "physical momentum" . The 
particle say for k > 1 has an energy e> (fc) which lies between and k 2 — k. We can 
also view the excited particle state as a combination of a particle and particle-hole pair, as 
follows. We write a schematic development for k > 1 

c\k) ~ ]T A t (p)A t (Jfe-p + 1)B\-1) , (49) 

i±*<p<fc 

The excitation energy of this complex is readily seen from Eqs. ( |39"| , [4~1~D to be EA(p) + EA{k — 
p + 1) + E B (-1), with ^ <p<k. The sum of these three terms reproduces the variation 
in e>\k\6) implied by the rapidity variable. The density fluctuation p q is then seen to be 
formally a two quasi particle-hole object: writing c(k) ~ B\—k), we have 

c\k)c{k-q)~ Y, A\p)A\k —p + l)B\-l)B\q - k) , (50) 

where it is understood here and elsewhere that when two momenta coincide (as they would 
in say [B^(— l)] 2 ), then these should be separated by the smallest non-zero wave number. 
The above scheme for the density fluctuation operator v{k)c(k — q) is indicated in Fig. 11. 
We may therefore regard the density fluctuation as being built up from a particular set of 
(non-interacting) pair states consisting of annihilating two particles at momenta k — q and 
1, and creating a pair with total momentum k + 1, distributed over all possible relative 
momenta with appropriate form factors. 

Symplectic Ensemble: Change of Variables 

We recall for the symplectic case, the Bethe energies e>(p) = p 2 + 2\p\ and e<(p) = 2p 2 + l 
(Eq. P^D). We now rewrite Eq. (|2B[) using the same variables as in the previous case 
(Eq. (f4"5|)). We find the result breaks up naturally into two pieces S a and Sb, with the 
second piece S& only arising for q > 2: 

S s (q, u) = S a (q, u) + % - 2) S b (q, u) , (51) 
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with 

a 2 r 1 r( 1 + u ) d? N 

S a (q,co) = \-\ du l Z 6(E a )^ , (52) 

lb Jmax(0,l-q/2) J2^u \J Z 1 — AlU U a 

and 

a 2 r° /•(!+«) dz Nh 

ID Jmax{- 1,1-9/2} -/O V-2 — 4u -L'fi 

where i£, iV and D are appropriately defined (see below). We write u — (l + k)/2 in S a and 
u = (I — l)/2 in Sb, in terms of which Eqs. fl5"2"|) and (|5*3] ) take a more natural form 

S a {q,u) = h dkn (k)[l-n (k + q)] I S(E a )-^ , (54) 



32^ " " /J V 2 ^) -2(1 + k) ~" "'A. 

and 

S h (q,u,) = ±- dl n Q (l)[l - n (l + q - 2)} / 5(E b )—^ . (55) 

32 J JO , /r2 i 0/1 _ 7^1 1A 



^ 2 + 2(1 - /) 

In Eq. (Q) we further introduce the rapidity variable 6 through z = 8(1 + A;)sinh 2 #, so 
as to eliminate the square root in the integrand. The result can be written compactly as 
follows 



n / v q 2 f „ [e>(k + g)-3} r*y/vm 

Sa(q,U) = IT /|*|<i 0*7 o n -l 7M2 / ^ 

2 4+ g >i [u; - 2g(l + fc)J 2 Jo 
x $(A(jfe, g) + 8(1 + fc) sinh 2 - w) . (56) 

We can perform explicitly the rapidity integrals and find the final result 

q / ,_q 2 f AU [e>(fc + g)-3] 

b a fru) dk [uJ - A (k, q )-(k-l) 2 + q 2 ] 2 

6(Lj-A(k,g)) 9(A(k,q) + (k-l) 2 -cj) 
x — - =- . (57) 

tJlu - A (A;, gtyw - A(k, q) + 8(1 + k) 
We next turn to the other piece for q > 2. With a = (q — 2) and introducing the rapidity 
variable <p through z = 8(1 — Z)sinh 2 (0), Eq. ( |5"5j ) is expressible in the form 



S b (q,uj) = %- I \i\<i dl 



[e>(/ + a)-3] /-inV2/(i-0 



o 



2 7^a>i - 2g(/ - l)] 2 
x 8(A(l,a)+ 8(1 -/) sinh 2 - w) . (58) 



Performing the rapidity integration we find 

q 2 f A1 [e>(/ + a) -3] 



S b (q,uj) = \ [ |i|<i dl 

4 J/+rv>1 



[a;- A(Z, «)-(/ + 1) 2 + g 2 ] 2 
x fl(n;-A(i,a)) g(A(/,cQ + (Z + l) 2 -^) 
- A(Z, a) y'cj - A(Z, a) + 8(1 - T) 



Bethe Quasi Particle-Hole Content: Symplectic Ensemble 

The energy conserving delta function in Eq. flSTp can be rewritten using an effective 
energy variable e" (k\9) = e < (k) — 8(1 + A;) sinh 2 6, which implies uj = e > (k+q)—e<(k\9). The 
energy {k\8) varies between the limits 2A; 2 + 1 (i.e. £ < (At)) and fc 2 + 2A; ( i.e. (/c) — (A; — l) 2 ). 
It is thus evident that we may interpret £< (k\9) as an effective hole, i.e a composite object. 
One possible way to decompose it is to write schematically for the bare annihilation operator 
a representation as a quasi-hole plus a quasi particle-hole pair: 

c(Jfe)~ ^ B\p — k - l)B\-p)A\l) . (60) 

k<P<^- 

The restriction on the range of p is such that we avoid double counting the pair and have a 
natural ordering of the two quasi-holes. The energy of the effective hole is then Es(—p) + 
Eb{p — k — 1) + Ea(1), with the constraint k < p < ^jl, which from Eq. (|39|j4l| ) reproduces 
the range required by the rapidity variation. Therefore, the operator p q is seen to be formally 
a two quasi particle-hole object, 

c*(k + q)c(k)~ ]T B\p-k-l)B\-p)A\l)A\k + q) . (61) 

k< P <^ 

This scheme for the density fluctuation is illustrated in Fig. 11. 

In the second piece of S s (Eq. ([58])), once again the energy conserving delta function can 
be rewritten introducing an effective energy variable for the hole £<(Z|0) = £<(/) — 8(1 — 
Z)sinh 2 0, which implies uj = £>(Z + a) — £<(Z|0). The effective hole energy £<(Z|</>) varies 
between £<(£) and £<(£) — (1 + I) 2 , corresponding to a predominantly left moving object, 



19 



and may be decomposed again into a hole and a particle-hole pair. Schematically we have 

c(Z)~ B\-p)A\-l)B\l+p-l) . 

1 -y-<p<i 



Using the quasi-energies Eqs. ( |39l , |4TD this complex has energy Eb(p)+Eb(1+p— l)+E / 



with the physical constraint ^y- < p < /, which reproduces the range implied by the variation 
of the rapidity. Owing to momentum conservation, we must regard the creation operator 
c*(l + q) as A\l + q — 2) times a particle- hole pair with energy zero and momentum 2, i.e. 
A'(1)B*(1). We may eliminate a 'zero pair' A\— 1)B<(1) and thus obtain the scheme for 
the density fluctuation operator p q , 

c\l + q)c(l)~ A\l + q-2)A\l)B\-p)B\l+p-l) . (62) 

The term Sb then evidently may be regarded as a two quasi particle-hole object, and is 
illustrated in Fig. 11. 

Summarizing, in process (a), we may regard the density fluctuation as being built up 
from a particular set of (non-interacting) pair states consisting of creating two particles at 
momenta 1 and k + q, and destroying a pair with total momentum k + 1, distributed over all 
possible relative momenta with appropriate form factors. Likewise, in process (b), we may 
regard the density fluctuation as being built up from a particular set of (non-interacting) pair 
states consisting of creating two particles at momenta 1 and I + q — 2, and destroying a pair 
with total momentum I — 1, distributed over all possible relative momenta with appropriate 
form factors. 



V. SPIN 1/2 HEISENBERG SYSTEMS 

In this section we will study the moments I n in two standard spin 1/2 Heisenberg an- 
tiferromagnetic systems, the 1/r 2 system and the Bethe chain. The Heisenberg spin chain 
model with a 1/r 2 interaction was introduced by Haldane 0] and Shastry ||. It is defined 
by the Hamiltonian 
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where <fi = T{ = 0, 1, L — 1 (L integer), and the spins are 1/2. In this case the natural 
operators one can use to introduce a dynamical structure function S(Q,u) are the Sf: we 
define a "charge operator" pi = (Sf + 1/2) and use (|) to write 

S d (Q,u) = 1 £ |H^o|0>| a *(w - + #o) , (64) 
with E v as the energy eigenvalues of the Hamiltonian and the lattice Fourier transform 

Sq = J2 S-e- lQr3 > (65) 
i=i 

where Q is the lattice momentum Q = (2ir/L) x integer. N is the number of spin deviations 
or the number of hard-core bosons. We set iV = dL, so that d — 1/2 for half filling. The 
Fermi momentum is then hp — vrrf, and we will scale Q = ndQ in order to compare with the 
continuum model results. 

We cannot calculate S d (Q,u>) directly, but we know some moments of this distribution, 
just as we did for the continuum model. As in the continuum, we restrict ourselves to 
n = 0, 1, and —1. 



A. Static Correlation Functions and the Zeroth Moment 



We begin with the zeroth moment, or the static structure factor 



I«(Q) = / S d (Q,u)duj. 
Jo 



(66) 



There is a remarkable theorem by Mehta and Mehta |E3| stating that the static correlation 
function is identical to that of the repulsive (fl = 4) continuum model in real space. This 



object was also calculated independently for half filling in Ref. [24], and the result (in Q 
space) is for the half filled case: 



(67) 
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in the scheme where we restrict \Q\ < it. 



The correlation function I is also available from Mehta and Mehta for arbitrary 
densities d < |. The density- density correlator can be written as 



(poPr) = d£r,0 + (1 - <M ^ [1 + D(r) 



(68) 



in a manner similar to Eq.flTBD. The function D is given |23] explicitly in the thermodynamic 
limit as 



DW = - 



sm(2irrd) 



2nrd 

Using the relation Sf = fa — ~, we find 



27rrd ^sint\ (2nrd) cos(27rr<i) — sin(27rr<i) 
1 o t_ T/ (27rrJ) 2 



(69) 



(S^) = - + d(l - d)(5 r , - 1) + (1 - S r>0 ) d 2 D(r). 



(70) 



Inverting the Fourier series, we have 



/ d (Q) = l + d £exp(iQr) D( 



r). 



(71) 



We may convert the sum over r to an integral, remembering that Q and Q + 2tt x integer 
are equivalent. We will work in the reduced zone scheme \Q\ < it, for which two cases may 
be distinguished, case (A) d < \ and case (B) \ < d <\. The correlations are given for 
Q > 0, and may be obtained for negative Q by using the evenness in Q of Iq- For case (A) 



d < \ we find 



J d (Q) = 1 + 9(And - Q)A(Q) 



A(Q) 



Q_ 

4nd 



Q_ 

8nd 



In 



Q 



2-ad 



(72) 
(73) 



and for case (B) i < d < ~ we find 



/ d (Q) = 1 + 



-27r + 47rd) A(2vr - Q). 



(74) 



It can be checked for d = 1/2 that Eq. ( [74] ) is identical with Eq. (p7|). For d < 1/4, the 
expression in Eq. ( f75f) is identical to the symplectic case Eq. (|13"D, apart from a scale factor 
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of d. Eq. (|74f) is in fact nothing but the Umklapp reduction of the continuation of Eq. (f73|), 
i.e. Q is allowed to extend up to Adn, and the part beyond it is declared to belong to Q — 2n, 
after subtracting unity from the structure function. 

The correlation function, in scaled variables, is given by (compare Eq. fl20|) ) 

"iM dQ 



D(r/d) 



i/d 2 



exp(iQirr) I^Qird) - 1 



(75) 



The theorem of Mehta and Mehta asserts the equality of the scaled correlation functions for 
/3 — 2, 4 for all integer r, i.e. for f = d x integer 



D{r^f/d) =C{r ->r/d). 



(76) 



B. Other Moments 

As opposed to the case when particles are in the continuum, the moment I\{Q) is not 
interaction independent for a system where the particles sit on a lattice: this is a well-known 
effect of the lattice systems with interaction ||. However, we can work out an expression 
for the 1/r 2 spin chain, using the usual definition as a double commutator. In the remaining 
part we will assume that d — 1/2 and write 



h(Q) = ([[Sq, H], S^q\) . 



(77) 



Calculating the commutator, we find 



= 7 ER«)[1 - cos(Q(r, -r,))](^) 



sin 0r 



[l-cos(Qr)](^>. 



Using Eq. (|B"TD, we rewrite this as 

h\Q) = %- E Mi-l*IA0£ 

Using the fact that, for < tt Q, 



sin 0r 



[1 — cos(/cr)] cos(fcr) 



(78) 



(79) 
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E 



sin (0r\ 



n 2 ( \k\ 

cos(fcr) = —(1 - 1/L 2 ) - 7r|Jfe| 1 - — 
3 \ 2ir 



(80) 



we find with [k] = (k — 2nm) [m integer \\[k}\ < ir] 

\[Q + k}\ + \[Q - k}\ - 2\k\ 



-^{[Q + A;] 2 + [g-A;] 2 -2A; 2 } 



(81) 



After turning the sum to an integral, we can integrate the expression and find, for \Q\ < tt, 

Jtt 2 



l 



\j3\\ 2 ( 1 m 2 hl ( 1 _m 
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7T 



TT 



(82) 



In the limit Q — > 0, we see that lf{Q) — > JQ 2 /4 but for larger Q there is substantial 
departure from the pure quadratic behavior of the continuum models Eq. (0). 

In the hydrodynamic limit the moment I_i(Q) can be obtained from the spin suscepti- 
bility Xspin, which is known explicitly ||: 



i d -AQ) 



L 



' d*E 
dM 2 



M=0. 



Xspin{Q > 0) 



(83) 



L 7T 2 J 

where M. = (N^ — iVj)/L is the magnetization and the operators count the number of 
up or down spins in the chain (notice that we have + — L/2). 

With these moments at hand, we can define characteristic frequencies for the excitation 
spectrum, in the same way as we did before. The important point is that all those frequencies 
will have the same value as Q — > 0: 



if(Q) Q^o JQ| 

i d i(Q) " ~\Hi{Q)^ 2 



(84) 



This means that in the hydrodynamic limit the spectrum is exhausted by the excitations 
with dispersion relation ui = Jtt\Q\/2 ||||. Although this might be taken as a common 
feature for continuum systems, this is not so in general, as we shall see below. 

We see that, by choosing appropriate energy scale (J — ► ^|), the 1/r 2 discrete spin model 
moments map exactly onto the continuum symplectic ones for Q — > 0. The dimensionless 
moments are identical in this limit, 
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iJ(g = 7rg/2)«il'(g = M), 

ftp 



(85) 



where the factor of 1/2 arises from the different normalization in Eqs. (§,0). It would be 
very interesting to pursue the calculation of more moments at q > for both models to 
check whether the discrete one shares the same unusual characteristics of the continuum 
(j3 = 4) model excitation spectrum. 

We now turn to the Bethe chain, which is defined by the Hamiltonian H = JnY^i SrSi+i- 



This model was studied in Ref. [25|, and does not show the saturation at Q — > 0. The 



moments I\ and J_i are known |p5 |, but not Iq\ 



jf(Q) = ~(21n2-l/2)(l-cosQ) 



(86) 



and 



i-AQ) 



n 2 J 



(87) 



The spinon spectrum of this system is known from Faddeev and Takhtajan's work 26 



^spiQ) — ( 7r ^/2) sin \Q\. Therefore, we could look for saturation by the spinon, and hence 
form the ratio 



lim 



1 



CO 



h B (Q) 



1.08706 . 



(88) 



It is therefore clear that the small Q behavior of S(Q, u) is not exhausted by the spinons in 
the Bethe chain, unlike in the 1/r 2 model, where it is. 



VI. CONCLUSIONS 

We have seen that the results of the calculation of Simons et al |Tl|-4T3| for the dynam- 
ical structure function has a representation in terms of the Bethe quasi-particle quasi-hole 
energies. The representation obtained in this work, has the character of two particle-hole 
pairs representing the bare density fluctuation. We should note, however, that this repre- 
sentation is far from unique, for one thing one may add an arbitrary number of "zero pairs" . 
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Also, for example, we could decompose our two holes and a particle as three holes and two 
particles, by e.g. forcing the momenta of two holes to coincide, and by suitably restricting 
the momenta. However, it is not possible to decompose the results into those of a single 
particle-hole pair; our representation in terms of two pairs appears to be minimal in some 
sense. The striking feature which underlines the dynamical structure factors presented in 
this paper is the truncation at very low orders of series like Eq. fl42|) . Therefore, the ex- 
cited states for system with (3 = 1,2, and 4 will always involve a small number of quasi 
particle- hole pairs. This may not be true for arbitrary values of the coupling constant. 

We have also shown that the discrete model shares the property of saturation of the 
Feynman sum rule by the lowest mode as q — > 0. The static structure function obtained 
by a direct calculation p4| is shown to be consistent with the older calculation of Mehta's 



Ref. |23| provided one interprets the weight outside the Brillouin Zone appropriately by 



Umklapping it. The first moment of the discrete model is obtained using the known result 
for the two-point static correlator, and shows interesting structure and departure from the 
first moment of the continuum model, and should provide a non trivial check on the structure 
function of the discrete model. 

We stress that the saturation of the structure function by the sound modes at small q is 
a general property characterizing this class of models. 
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FIGURES 

1. Characteristic frequencies of S(q,u). The solid line corresponds to the Feynman 
spectrum u)p(q) and the dashed line corresponds to the Hydrodynamical spectrum u-n(q): 
(a) (3 = 2; (b) (3 = 1; (c) (3 = A. 

2. Region where S u (q, uo) ^ {(3 = 2). The equations for the boundaries are indicated. 

3. Tridimensional plot of S(q,u) {(3 = 2) for the unitary case. The vertical axis has a 
linear scale. 

4. Moments of S(q,u) for (3 = 2: (a) I x {q)] (b) I {q); (c) 

5. Region where S°(q,u>) ^ (/3 = 1). The equations for the boundaries are indicated. 

6. Tridimensional plot of S(q, u) for the orthogonal case (fl = 1). The delta function of 
Eq. ( f24"|) has been regularized by a Gaussian. The vertical axis has a logarithmic scale. 

7. Moments of for (3=1: (a) Ji(g); (b) I (q); (c) J_i(g). 

8. Region where S s (q,u) ^ (/3 = 4). The equations for the boundaries are indicated. 

9. Tridimensional plot of S(q,u>) for the symplectic case (/3 = 4). The delta function of 
Eq. ( ^6|) has been regularized by a Gaussian. The vertical axis has a logarithmic scale. 

10. Moments of S(q,u) for (3 = A: (a) ^(g); (b) J (g); (c) I-x{q). 

11. The two quasi-particle quasi- hole pair scheme for the orthogonal and symplectic 
cases. 'X' denotes particles and 'O' denotes holes, and the solid line indicates the range 
— 1 < k < 1, i.e. the Fermi sea (kp = 1). The symplectic ensemble has two pieces: case (a) 
corresponding to S a and case (b) corresponding to Sb- 
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